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ABSTRACT 


Impact-induced devolatilization of porous serpentine was investigated using two independent 
experimental methods, the gas recovery method and the solid recovery method, each yielding 
nearly identical results. For shock pressures near incipient devolatilization CPinitial = 5.0 GPa), the 
hydrogen isotopic composition of the evolved HjO is very close to that of the starting material. For 
shock pressures at which up to 12% impact- induced devolatilization occurs (up to Pjnjtial = 7.7 
GPa), the bulk evolved gas is significantly lower in deuterium than the starting material. There is 
also significant reduction of H 2 O to Hj in gases recovered at these higher shock pressures, 
probably caused by reaction of evolved H 2 O with the metal gas recovery fixture. Gaseous H 2 O-H 2 
isotopic fractionation suggests high temperature (1 100-1200 K) isotopic equilibrium between the 
gaseous species, indicating initiation of devolatilization at sites of greater than average energy 
deposition (shear bands). Bulk gas-residual solid isotopic fractionations indicate nonequilibrium, 
kinetic control of gas-solid isotopic ratios. Impact-induced hydrogen isotopic fractionation of 
hydrous silicates during accretion can strongly affect the long-term planetary isotopic ratios of 
planetary bodies, leaving the interiors enriched in deuterium. Depending on the model used for 
extrapolation of the isotopic fractionation to devolatilization fractions greater than those investigated 
experimentally, planetary interior enrichments of 4 to 30 per cent in D/H relative to that of the 
incident material can result from this process. 
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INTRODUCTION 


Impact-induced devolatilization of volatile-bearing rocks and minerals is an important process 
influencing planetary accretion, early atmosphere formation, and chemical weathering of planetary 
surfaces [1,2,3,4,5J. Previous experimental investigations of this process have focussed on the 
shock pressure (impact velocity) dependence of the extent of impact-induced devolatilization in 
hydrous minerals, carbonates, and meteorites [6,7,8,9,10] and on the shock enhancement of 
chemical reactivity of silicates [1 1,12]. In this study, we perform a direct comparison of the two 
experimental techniques for measuring impact-induced devolatilization, namely, the gas-recovery 
and solid-recovery techniques. In addition, we have examined the hydrogen isotopic composition 
of the gaseous products of impact-induced devolatilization of serpentine. 

Early experimental investigations of impact-induced devolatilization employed gas recovery 
techniques, in which the gas driven from a sample during impact was recovered directly and 
analyzed gravimetrically [6] or manometrically and isotopically [7]. The purpose of these early 
studies was to establish the existence of the impact-induced devolatilization process. Accordingly, 
shock pressures were investigated only up to the point of incipient devolatilization, the shock 
pressure at which the material just begins to release volatiles upon impact. Subsequent studies that 
more fully delineated the shock pressure dependence of the extent of impact-induced 
devolatilization utilized solid recovery methods, in which solid material was recovered after the 
impact and analyzed for its post-shock volatile content [8,9,10]. Differences in volatile content 
between unshocked material and recovered, shocked material were attributed to the impact process. 
A potential difficulty with the solid recovery method is the problem of the handling of the materials 
between the shock event and the volatile analysis. The shocked assemblies must be machined open 
to recover the sample, and given the extremely reactive nature of these materials, the question can 
be raised whether even brief exposure to the ambient atmosphere compromises the samples. 

This study was performed 1) to directly compare impact-induced devolatilization results using 
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gas recovery and solid recovery techniques, and 2) to examine the chemical and isotopic 
composition of the gases evolved during impact into a hydrous sample. Direct comparison of solid- 
and gas-recovery techniques has not been made previously, although the higher pressure solid 
recovery data are consistent with the incipient devolatilization pressures of the earlier gas-recovery 
work. In the only study to date of the isotopic composition of shock-evolved gases, CO 2 evolved 
by shocking calcite just to the point of incipient devolatilization had the same C- and O-isotopic 
composition as the starting material [7]. On the other hand, it has been suggested that vapor 
fractionation of impact-produced melt could provide an explanation for the relatively low O 
ratios in tektites [13]. No previous examination of the isotopic composition of gases driven by 
impact from a hydrous mineral has been performed. 

EXPERIMENTAL 

The impact experiments were performed on a nearly pure magnesium end-member, 
polyciystalline lizardite serpentine collected near Globe, Arizona (Caltech Mineral Collection, 
Specimen Number 3221). The bulk chemical analysis is shown in Table 1. The density of the 
mineral is 2.50 ± 0.02 g/cm^; 8 Dj M q w is -77 per mil. 

Solid-recovery impact-induced devolatilization experiments were performed on the 20 mm 
powder gun at the California Institute of Technology using vented, stainless steel sample 
assemblies and thermogravimetric analysis of the recovered material as described previously 
[6,10]. Sample powders pressed to approximately 80 to 85 % of the bulk mineral density were 
studied. Porous materials were studied because the high temperatures accompanying shock 
compression of porous materials lead to increased devolatilization at a given shock pressure relative 
to nonporous materials. Thermogravimetric analysis (TGA) of shocked and unshocked samples 
was performed over the temperature range 300 to 1 100°C using a Mettler Thermoanalyzer 2000C. 
Sample size was approximately 20 mg. Analysis was performed in a nitrogen atmosphere. The 
extent of impact-induced devolatilization was determined by comparison of TGA analyses of 
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shocked and unshocked material. Table 2 lists the results of the solid recovery experiments. 

Gas-recovery experiments were performed using the assembly sketched in Figure 1. This 
assembly is a modification of that employed by Boslough et al. [7]. The sample powder was 
pressed into a wafer, inserted into the assembly, and held in place by the threaded containment 
plug. Sample densities were approximately 80 to 85 % of the bulk mineral density. Sample 
dimensions were about 0.5 mm in diameter by 0.5 mm thick. Sample masses were between 14 and 
21 mg. The impacted surface-to-sample distance was approximately 0.6 cm. After loading the 
sample, the chamber was assembled, evacuated to less than 10 pm Hg pressure using a diffusion 
pump, and heated at 75°C overnight Prior to the impact experiment, the assembly was again 
evacuated, the valve tighdy closed, and the assembly mounted on the gun barrel. After the impact, 
the assembly was connected to a vacuum system, and the trapped gases were transferred through 
the valve section to the vacuum system for manometric and isotopic analysis. The integrity of the 
sample assembly was maintained for impact velocities up to 1.22 km/s. Impact velocities greater 
than 1.3 km/s resulted in rupture of the chamber. A blank experiment was run in which no sample 
was placed into the assembly, and the recovered gases were analyzed in the same manner as the 
sample-bearing experiments. In general, in the gas recovery experiments the solid could not be 
recovered; it was spread as a fine dust over the interior of the assembly. In one shot, however, 
enough solid material was recovered to perform TGA (shot 867). The independent results for this 
shot are listed in Tables 2 and 3. For this shot, the agreement between the two techniques is 
excellent 

Manometric analysis of the recovered gases was performed using standard techniques [14]. 
The first separation was into noncondensible and condensible fractions at liquid nitrogen 
temperature. The noncondensible portion was further separated into a fraction that condensed onto 
molecular seive at liquid nitrogen temperature and the portion that did not condense. Subsequent 
mass spectrometric analysis of this 'doubly' noncondensible material established its identity as H 2 , 
and its D/H ratio was determined. CO 2 was released from the first condensible portion with an 
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organic solvent-Diy Ice bath. The final H 2 O portion was reduced to H 2 by passing it through 
uranium at 700°C and then analyzed mass spectrometrically. The results of the manometric and 
isotopic analyses for the blank- and sample-derived gases are listed in Table 3. Figure 2 illustrates 
the amounts of each gas fraction recovered for each of the gas-recovery experiments. The 
background levels of non-condensible gases, H 2 , and H 2 O recovered are significantly lower than 
amounts recovered from sample-containing shots, indicating that these recovered gases originate 
from the sample. Recovered COj is on the order of 1 |imole or less in all of the shots. The blank, 
impacted at 1.11 km/s, yielded 0.4 pmoles of H 2 and 1.0 |imole of H 2 O, amounts insufficient for 
accurate mass spectrometric determination of D/H ratios. These amounts were subtracted from the 
total H 2 plus H 2 O recovered in determining the fraction of initial H 2 O volatilized during the 
impact 

Shock pressures were calculated using a one-dimensional impedance match solution [15]. We 
report both the initial shock pressure and the peak reverberated shock pressure in Tables 2 and 3. 
Equation of state parameters for stainless steel were taken from [16]; for serpentine we used initial 
(nonporous) density p 0 = 2. 50 g/cm^, C 0 = 5.30 km/s, and s = 0.78 [17], where shock velocity 
U s and particle velocity Up are related by U s = C 0 + s Up. The Hugoniot of porous serpentine was 
calculated using the relation 


p H,p = p H,d d - YKPh'Po) ' 1 3 /2 > / d ‘ WPtfPoo) - (D 

[18]. In equation (1) Pn,p * s Hugoniot pressure of the porous material, Pfj.d ^ Hugoniot 
pressure (at ppj) of the nonporous material, y is the Griineisen parameter, pjj is the Hugoniot 
density, p c is the initial density of the nonporous material, and p^ is the initial porous material 
density. The Griineisen parameter is calculated using y = y 0 ( p 0 /Pj|) n Yo = ^ n = 1*0 

[17]. Initial densities, impact velocities, and calculated initial and peak shock pressures are listed in 
Tables 2 and 3. Experimentally determined errors in p 0 , p^, and projectile velocity Vp were 
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assumed to be independent and were propagated through the impedance match relations and 
equation (1) to yield uncertainties in initial and peak shock pressures [19]. 

RESULTS AND DISCUSSION 

For the higher shock pressures studied, both H 2 and H 2 O are found in the recovered gases and 
the D/H ratios of both the H 2 and H 2 O produced are significantly lower than that of the starting 
material (compare shot 867 with shots 911 and 915,Table 3). We postulate that the H 2 in the 
recovered gas is a result of oxidation of iron metal in the experimental assembly according to the 
reaction 


Fe + H 2 0 (gas) = FeO + H 2 ( gas ) 


( 2 ) 


Under standard state thermodynamic conditions the products of this reaction are favored over 
reactants up to temperatures in excess of 1800 K [20], and for conditions in which Ph 20 > ^H2> 
as is the case for the measured gas compositions of this study, this tendency will be enhanced. 
Reduction of H 2 O to H 2 during pyrolysis of metal-containing carbonaceous chondrites has been 
observed [21]. In addition, measureable amounts of both H 2 and H 2 O have been observed in mass 
spectrometric analyses of the spreading evaporation cloud of laser-shocked carbonaceous 
chondrites [22]. 

We compute the total H 2 O driven from serpentine by the shock as the sum of the H 2 O plus H 2 
collected, less the blank (shot 913). As indicated in Table 3, the actual H 2 content of the 
noncondensible recovered gas fraction was not determined for each experimental shot On the basis • 
of the results for shots 911 and 915, we assume that approximately 1/2 of the total noncondensible 
gases recovered is H 2 . Accordingly, for shots 875 and 890, the total H 2 plus H 2 O recovered is 
calculated as the sum of the H 2 O recovered and 1/2 of the noncondensible gases recovered. Figure 
3 compares the impact-induced water loss for porous serpentine as a function of peak shock 
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pressure determined in this study by the gas- and the solid-recovery methods. In the pressure 
region of overlap, the agreement is excellent. We conclude that the solid and gas recovery methods 
provide independent measures of the same phenomena and that previous impact-induced 
devolatilization results based on solid recovery methods are valid. 

Figure 3 also shows the previous results for impact-induced devolatilization of nonporous 
(density 2.5 g/cm^) serpentine [8]. The shock pressures for this previous work have been 
recalculated using the more recent serpentine equation of state data. The porous material shows 
significantly enhanced devolatilization over the nonporous material at the same shock pressure, 
consistent with expectations based on the higher shock temperature and shock entropy content of 
porous material [23], The empirical best-fit polynomial curves for each data set have different 
curvatures. We interpret this difference as an artifact resulting from different shock-pressure 
distributions of the data in each set. That is, there are more low-pressure, low devolatilization 
fraction data for the porous material and more high-pressure, high devolatilization data for 
nonporous serpentine. The shock entropy model for impact-induced devolatilization indicates that 
the two curves should be of similar shape. 

The hydrogen isotopic composition (6D = {[(D/H) sample /(D/H) SM0W ] - 1} x 1000, where 
SMOW refers to Standard Mean Ocean Water) of the gaseous H 2 O recovered from shot 867 
^initial - 5-0 GPa) is nearly identical to that of the starting material, -72 compared to -77 per mil, 
respectively (Tables 1 and 3). However, for shots 911 and 915 OPjnitial = ^-4 and 7.7 GPa, 
respectively) the recovered H 2 O has a significantly lower deuterium concentration than the starting 
serpentine, -93 and -98 per mil, respectively. The isotopic fractionation factor for hydrogen 
between the recovered H 2 O and H 2 , a H20-H2 = (1000 + ^D^O^QOOO + SDj^) = 
(D/H)jj 2 Q/(D/H)pj 2 , and the per mil difference, Apj2Q-H2 = 10^ 1° a H20-H2> 816 listed in Table 
4. For these two shots, the Ajj20-H2 va l ues are 131 and 169 per mil, respectively. Hydrogen 
isotopic fractionation between gaseous H 2 O and H 2 has been measured to 1000 K [24] and 
calculated theoretically [25]. For the experiments reported here, initial and peak shock temperatures 
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were calculated using the procedures outlined by Tyburczy and Ahrens [23], the equation of state 
parameters for serpentine listed above, and these additional EOS parameters; zero-pressure 
adiabatic bulk modulus K os = 64 GPa, pressure derivative of the adiabatic bulk modulus K' = 2.7, 
and specific heat capacity at constant volume C v = 1.3 J/(g-K) [17].The calculated initial and peak 
shock temperatures lie between 600 K and 770 K (Table 4). Figure 4 is a plot of H 2 O-H 2 
deuterium fractionation ctH 20 -H 2 versus temperature, showing both the equilibrium curve and the 
experimental results of this study. The values for the impact-derived gases lie at significantly lower 
temperatures than the equilibrium curve. If one assumes that the impact-produced H 2 O and H 2 
were in isotopic equilibrium, the measured values of a H20-H2 indite equilibration temperatures 
of approximately 1100-1200 K. This result suggests that shock-induced devolatilization is initiated 
at localized sites of thermal energy deposition, such as 'shear instability' regions [26,27,28]. 

We calculate the bulk 8 D of the released gases from the 8 D and the molar amount of each gas 
species recovered. Then, using the initial 8 D of the serpentine and the fraction of H 2 O released 
during the shock, we obtain the bulk 8 D of the residual solid and the fractionation factor 
a gas-residue 311(1 ^ difference A gas _ rcsidue (Table 4). Values obtained for a gas . res idue 

are 1.005, 0.941, and 0.904 for shots 867, 91 1, and 915, respectively. Values of otpEQ-residue 
calculated using the 8 D composition of the H 2 O recovered lie between 1 .005 and 0.966 (Table 4). 
Figure 5 shows these results plotted versus shock temperature along with ®H20-mineral ^°r 
serpentine, muscovite, and rhyolitic melt for comparison [29,30,31,32]. Electron microscopic 
images of shocked serpentine reveal vesicle-containing glassy regions, suggesting that the volatiles 
are evolved from shock-melted regions of the sample [ 8 ]. The trend of <*H20-residue is 
inconsistent with the temperature trends of equilibrium I^O-mineral or I^O-melt D-H 
fractionation. The trend of 04 ^^ g as - re sid U e versus shock temperature (Table 4) would be even 
more discordant with the equilibrium curves. We conclude that the observed gas-residue isotopic 
fractionations are the result of kinetic processes occurring during the impact event and do not reflect 
equilibrium fractionation. 
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Thus, for shock pressures CPinitial = 5-0 GPa) near the point of incipient devolatilization, there 
is little or no reduction of H 2 O to H 2 in the shock-evolved gases, and the hydrogen isotopic 
composition of the gas is very close to that of the starting material. For higher shock pressures, at 
which more impact-induced devolatilization occurs (Pjnitial ^ 6.4 GPa), the bulk evolved gas is 
significantly lower in deuterium than the starting material and there is significant reduction of HjO 
to H 2 , probably caused by reaction of evolved H 2 O with the container. The lack of isotopic 
fractionation at lower shock pressures is consistent with the results of Boslough et al. [7] for 
calcite. They found only minor variations in 8^C and 8^0 between the starting material and the 
CO 2 evolved for calcite shocked to incipient devolatilization pressures. 

The influence of impact-induced hydrogen isotopic fractionation of hydrous silicates during 
accretion on planetary isotopic compositions depends on the mass fraction of the planet accreted in 
an impact velocity regime in which partial devolatilization of incident and impacted material occurs. 
For very high impact velocities, complete devolatilization occurs and no isotopic fractionation 
would take place. In order to model the effects of impact-induced hydrogen isotopic fractionation 
on planetary accretion, information on the fractionation over the entire devolatilization interval is 
required. Because of the impact-velocity limitations of our gas-recovery chambers, the present 
results extend only to devolatilization fractions of 12%. Figure 6 shows the experimental bulk 
gas-residual solid fractionation factor a bu y c g as -residue as a function of devolatilized fraction Xj} 
results of this study along with model extrapolations to higher Xq. Lacking a firm basis for 
extrapolation, we examine three possible models. Model 1 is a simple linear extrapolation of the 
experimental results; oc^]^ gas-residue reac h es a minimum value of 0.33 at Xj) = 1.0. In model 2 
we assume that gas-residue decreases linearly to 12 % devolatilization, then remains constant * 
at a value of 0.9. In model 3 we assume a simple parabolic function through the data, reaching 

a bulk gas-residue = at 

To model the effects of impact-induced hydrogen isotopic fractionation on planetary accretion, 
we assume a very simplified homogeneous accretion model. We assume that the incident material 
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and the planetary surface both are composed of serpentine (equation of state discussed above), and 
that the growing planet maintains a constant density equal to its final average density throughout the 
growth period. Under these conditions the minimum impact velocity Vp is a linear function of 
radius r 

V p = r-(8 7tGp avg /3) 1/2 (3) 

where G is the gravitational constant and p avg is the (present-day) average density of the planet 
[10]. Figure 7 is a plot of impact-induced devolatilization as a function of fractional radius for the 
Earth, Venus, Mars, and the Moon, using the parameterized devolatilization versus initial shock 
pressure expression for nonporous serpentine (dashed curve in Figure 3). For the Earth and Venus, 
impact-induced devolatilization of incident planetesimals would commence when each planet 
reached about 12% of its final radius; complete devolatilization would occur by the time the 
growing planet reached about 30% of its final radius. For Mars, incipient and complete 
devolatilization would occur at about 28% and 70%, respectively, of the final radius. For a planet 
of the Moon's size and density, impact induced devolatilization begins at 60% of final radius, and 
at its final size, incident material would be 60% devolatilized on impact The results presented in 
Figure 7 differ slightly from previous calculations [10] because of updates in the serpentine 
equation of state parameters. 

We also assume that the volatile materials not released on impact remain associated with the 
solid planet and are subsequently buried. The volatiles released on impact either remain in the 
planetary atmosphere or recondense on the planetary surface and are remobilized during subsequent 
impacts. These simplifying assumptions allow us to calculate the total surface volatile budget or 
cumulative fraction devolatilized for each of these planetary bodies (Figure 8). For the Earth and 
Venus, nearly 99% of the incident volatiles remain at the planetary surface at the end of accretion, 
for Mars 86% remain at the surface, and for the Moon, 27% remain at the surface. In each case, the 
remainder is retained in the planetary interior. 

For a given impact, the D-H fractionation relative to the incident material is calculated as 
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follows. Mass balance requires that 

(D/H)j = X D -(D/H) gas + (1 - X D )-(D/H) solid (4) 

where (D/H)j is the D/H ratio of the incident material and Xj) is the fraction devolatilized 
Substituting (D/H) gas = a bu]k gas-residue'^^solid ^ rearranging yields the solid-initial 
fractionation 

^solid-initial = ^^^solid^^^^i = + ^D*( a bulk gas-residue " ^ ($) 

Figure 9 shows the bulk gas and solid fractionations for Mars as a function of planetary radius, 

illustrating the effects of the different a bu lk gas-residue'^D m odels (Figure 6). The linear model, 
the most extreme of the models, yields a maximum (D/H) sobd nearly 3 times (D/H)j occurring at 
the radius at which complete devolatilization occurs. The minimum ^solid-initial * s approximately 
0.8, occurring at r/R = 0.5 when Xq = 0.62. The more conservative models 2 and 3 reach 
maximum a so lid-initial va ^ ues °f about 1.1 and minimum values near 0.9. 

The cumulative effects of these fractionations are shown in Figure 10. Assuming that all the 
H 2 O devolatilized on impact accumulates in the atmosphere or on the surface, and that the 
remainder remains within the solid planet, cumulative bulk solid-incident fractionations at the end 
of accretion (r/R = 1.0) are 1.296, 1.043, and 1.061 for models 1, 2, and 3, respectively. 
Cumulative bulk gas-incident fractionations are 0.96, 0.99, and 0.99, respectively. Thus, 
depending on the exact nature of the gas-solid fractionation versus devolatilization relation, 
impact-induced hydrogen isotopic fractionation can lead to significant solid planet-atmosphere 
isotopic differences. 

In Figures 1 1 and 12 we compare the effects of impact-induced isotopic fractionation using 
model 3 (the conservative, parabolic model) on the development of the Earth, Venus, Mars, and the 
Moon. For the Earth, Venus, and Mars, each planet reaches sufficiently great size that complete 
devolatilization occurs during accretion (Figure 7). In each of these cases, the cumulative bulk solid 
D/H ratio reaches the same final value of 1.06, enriched in D relative to the incident material (Figure 
12), because beyond the critical mass (radius) for complete devolatilization no additional volatile 
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material is added to the planetary interior. The cumulative bulk atmospheric D/H ratio varies 
continuously over accretion. For the Earth and Venus it asymptotically returns to that of the incident 
material. The atmospheric D/H ratio of Mars and a Moon-sized object would be 96% and 89%, 
respectively, of the that of incident material. 

Long term persistence of the impact-induced isotopic fractionation signature would further 
require isolation of the gas and solid isotopic reservoirs, either by rapid burial of the solid shocked 
material or by rapid planetary escape of the evolved gases. In the event of rapid removal 
(catastrophic or otherwise) of the earliest planetary atmosphere, the bulk planetary D/H ratio would 
be the cumulative solid D/H ratio and would therefore be greater than that of the incident material. 

The cumulative D/H fractionations resulting from isotopic fractionation upon impact-induced 
devolatilization during accretion are not sufficient to explain the observed ~6- and ~ 10-fold 
enrichments in the D/H ratios of the atmospheres of Mars and Venus relative to the Earth [33,34]. 
However, this previously unrecognized process can cause different reservoirs within an accreting 
planet to have D/H ratios significantly different than those of the material from which the planet 
accumulated. The existence of these distinct reservoirs at the earliest stages of accretion will 
strongly influence the isotopic evolution of the body. 

Note that the experiments reported here and the above discussion deal with fractionation of 
hydrogen from hydrous silicates. Hydrous silicates in meteorites have 8D values around -100 per 
mil. An extremely wide range of 5D values (-500 to +10,000 per mil) is found for the 
acid-insoluble organic material of chondritic meteorites [21,35]. The shock behavior of these 
materials is not known, but their occurrence in accreting materials would significantly affect D/H 
ratios in a growing planet. 


CONCLUSIONS 

Two independent methods for investigating impact-induced devolatilization, the gas recovery 
method and the solid recovery method, yield nearly identical results for porous serpentine. For 
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impacts into porous serpentine near the point of incipient devolatilization (Pineal = 5.0 GPa), the 
evolved hydrogen-containing gases consist nearly entirely of H 2 O with a hydrogen isotopic 
composition very close to that of the starting material. For shock pressures at which up to 12% 
impact-induced devolatilization occurs (P^^ = 67.7 GPa), the bulk evolved gas is significantly 
lower in deuterium than the starting material and there is significant reduction of H 2 O to H 2 , 
probably caused by reaction of evolved H 2 O with the container. Gaseous H 2 O-H 2 isotopic 
fractionations suggest high temperature (ca. 1100-1200K) isotopic equilibrium between the 
gaseous species, but bulk gas-residual solid fractionations indicate nonequilibrium, kinetic control 
of gas-solid isotopic fractionation. Impact-induced hydrogen isotopic fractionation of hydrous 
silicates during accretion can affect the long-term planetary D/H ratios of smaller planetary bodies. 
D/H ratios of planetary interiors may be enriched by factors of 5 to 30 per cent over that of the 
incident material, depending on the exact nature of the isotopic fractionation-impact devolatilization 
relation. 
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Table 2. Solid Recovery Data 


Shot Number Initial Sample 



Density 


g/cm3 

867 

2.09 

±0.04 

873 

2.05 

±0.04 

870 

2.03 

±0.04 

869 

2.05 

±0.04 


Projectile 

Initial Shock 

Velocity 

Pressure 

km/s 

GPa 
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±0.01 
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1.27 
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±2.3 


Peak Shock 

Impact-induced 

Pressure 

H20 Loss 

GPa 

% 

17.5 

3.9 

±0.2 


27.7 

10.6 

±0.3 


36.4 

38.3 

±2.7 


43.8 

89.0 

±0.3 
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Shot Number Initial Sample Sample Mass Projectile Initial Shock Peak Shock Recovered Recovered Recovered Recovered Recovered Recovered Total Recovered 
Density Velocity Pressure Pressure Non-condensible H2 H2 C02 H20 H20 H2 ♦ M20 

Gases as % H20 in 

g/cm3 mg km/s GPa GRa micromoles micromoles dD, per mil micromoles micromoles dD, per mil Inital Sample 



- Recovered II 2 calculated as 1/2 total non-condensibles 



Table 4, Derived isoptop ic parameters 
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FIGURE CAPTIONS 


Figure 1. Schematic diagram of the gas recovery assembly. Valve, Be-Cu alloy; all other parts, 
stainless steel 304. 

Figure 2. Bulk chemical composition of recovered gases evolved from serpentine on impact. Shot 
913 was an empty, background experiment impacted at 1.1 1 km/s. Other shots listed in order of 
increasing impact velocity. 

Figure 3. Impact-induced water loss as a function of initial shock pressure for 20% porous 
serpentine (this study) and nonporous serpentine [8]. Filled squares, solid recovery method, 
porous serpentine (Table 2); open squares, gas recovery method, porous serpentine (Table 3); 
crosses, solid recovery method, nonporous serpentine. Solid line, empirical fit to data for porous 
serpentine; dashed line, empirical fit to data for noporous serpentine. 

Figure 4. H 2 O-H 2 D/H fractionation factor versus temperature. Curve is theoretical result of [25]. 
Data are results for shots 911 and 915 plotted against calculated shock temperature (see text ). 

Figure 5. I^O-residual solid deuterium fractionation versus shock temperature (Table 4), open 
squares. Solid curves, I^O-mineral deeuterium fractionation [29,30,31,32]. 

Figure 6. Impact-induced bulk gas-residual solid fractionation gas-residue 38 a f unc ti° n of 
devolatilization fraction Xp. Curves are assumed models for extrapolation (see text). 

Figure 7. Impact-induced devolatilization of nonporous serpentine as a function of fractional 
planetary radius. 

Figure 8. Cumulative impact-induced devolatilization of nonporous serpentine as a function of 
fractional planetary radius. 

Figure 9. Impact-induced D/H fractionation during accretion for Mars. Numbers refer to models of 
Figure 6. 

Figure 10. Cumulative impact-induced D/H fractionation during accretion for Mars. Numbers refer 
to models of Figure 6. 

Figure 11. Impact-induced D/H fractionation during accretion for Earth, Mars. Venus, and the 
Moon, using model 3 for bulk gas-residual solid D/H fractionation (Figure 6). 



Figure 12. Cumulative impact-induced D/H fractionation during accretion for Earth, Mars, Venus, 
and the Moon, using model 3 for bulk gas-residual solid D/H fractionation (Figure 6). 
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